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We consider an unambiguous identification of an unknown coherent state witli one of two unknown 
coiierent reference states. Specifically, we consider two modes of an electromagnetic field prepared 
in unknown coherent states |qi) and \oi2), respectively. The third mode is prepared either in the 
state or in the state \a2)- The task is to identify (unambiguously) which of the two modes 
are in the same state. We present a scheme consisting of three beamsplitters capable to perform 
this task. Although we don't prove the optimality, we show that the performance of the proposed 
setup is better than the generalization of the optimal measurement known for a finite-dimensional 
case. We show that a single beamsplitter is capable to perform an unambiguous quantum state 
comparison for coherent states optimally. Finally we propose an experimental setup consisting of 
2N — 1 beamsplitters for unambiguous identification among unknown coherent states. This setup 
can be considered as a search in a quantum database. The elements of the database are unknown 
coherent states encoded in different modes of an electromagnetic field. The task is to specify the 
two modes that are excited in the same, though unknown, coherent state. 



PACS numbers: 03.67.Dd, 03.65.Yz, 03.67.Mn, 02.50.Ga 



I. INTRODUCTION 

The ability to discriminate quantum states plays an 
important role in quantum information processing. Be- 
cause of the quantum interference two (non-orthogonal) 
quantum states cannot be distinguished perfectly pro- 
viding the number of copies of these states is limited. 
The topic of quantum state discrimination was firmly es- 
tablished in 1970s by pioneering work of Helstrom 
who considered a minimum error discrimination of two 
known quantum states. In this case the state identifi- 
cation is probabilistic. Another equally significant ap- 
proach is the unambiguous discrimination of quantum 
states, originally formulated and analyzed by Ivanovic, 
Dicks and Peres 0, 0, [1] in 1987. In contrast to the min- 
imum error discrimination approach, the unambiguous 
state identification is deterministic, i.e. no erroneous con- 
clusions are permitted. But in addition an inconclusive 
result is allowed corresponding to situations in which the 
state identification fails. The solution for unambiguous 
discrimination of two known pure states appearing with 
arbitrary prior probabilities (further denoted as r]i,r]2) 
was obtained by Jaeger and Shimony Q . Subsequent re- 
search was mainly focused on unambiguous discrimina- 
tion among several known pure states and unambiguous 
discrimination of two mixed states, which is still an open 
problem. The physical implementation of the optimal 
unambiguous discriminator device working for arbitrary 
coherent states was proposed by K. Banazsek Q. 

All results mentioned above are heavily based on a 
prior classical knowledge we have about those quantum 
states that are to be discriminated. S. Barnett et al. 
0] have studied this problem in detail and in addition 



they have addressed the following intriguing question: Is 
it possible to say anything unambiguously whether pure 
quantum states of a pair of identical quantum systems (fi- 
nite dimensional) are equal or not? Here no prior knowl- 
edge about the states is assumed. This problem is called 
quantum state comparison (for extension to more sys- 
tems see Ref. Q). It turned out that symmetry with 
respect to the exchange of the subsystems enables one to 
reveal unambiguously the difference between the states 
of the subsystems, however their equality cannot be de- 
termined unambiguously. E. Andersson et al. Q investi- 
gated this problem also for coherent states and proposed 
a simple setup (described in Sec. II), consisting of one 
beam splitter and a photodetector, capable to perform 
this task. 

Such results stimulated J. Bergou and M. Hillery [l3| 
to reduce the prior knowledge in unambiguous discrimi- 
nation of quantum states. These authors formulated the 
following problem: Imagine we are given two qubits B 
and C each of them in an unknown pure state. At the 
same time we are given also a third qubit A, which is 
guaranteed to be either in the state of the first or the 
second qubit. The task is to determine unambiguously 
with which of the two qubits the state of the third qubit 
matches. In such modification of the original problem 
the whole information is conveyed by states of quantum 
systems. Given the fact that we have just a single copy 
of each of the two states even the optimal quantum me- 
chanical measurement would allow us to determine those 
states with a rather small fidelity (for more details see 
Refs. P, [ll|, Therefore, instead of trying to esti- 

mate reference states we directly use them in their "quan- 
tum" form. The control of the behavior of the device by 
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the quantum program register is the key feature of pro- 
grammable quantum devices [H, 0, (for a review see 
Ref. [lB|)- In these devices one subsystem serves as the 
data register and the second quantum subsystem as the 
program register that carries instructions about transfor- 
mations the machine has to perform on the data register. 

The problem of unambiguous discrimination of un- 
known states, called also the unambiguous quantum state 
identification problem (UI), has been investigated over 
last few years by many authors. Bergou et al. [l8| exam- 
ined the situation with more copies of the third qubit A. 
If qudits [d dimensional quantum systems) instead of the 
qubits arc used and more copies of qudits A and B arc 
provided then the analytical solution for equal prior prob- 
abilities was obtained by A. Hayashi et al. in Ref. p^ . 
C. Zhang and M. Ying [20| investigated the unambiguous 
identification among n unknown qudit states and they 
provided necessary and sufficient criterion characterizing 
all possible programmable discriminators performing this 
task. 

One of the aim of this paper is to illustrate that a 
prior knowledge of a subset of states S uniformly en- 
tering subsystems A and B can significantly affect the 
optimal UI measurement and the performance of it with 
respect to the universal UI measurement. We illustrate 
this on two examples: The first example deals with the so 
called equatorial qubits described in Sec. II. B. The sec- 
ond example deals with coherent states examined in Sec. 
IV. In Sec. Ill we present how an "intuitive" universal 
UI measurement of systems of an arbitrary dimension 
can be constructed. Further we will compare different 
types of UI devices with the optimal universal one found 
in Ref. [l^. Although the UI measurement of coherent 
states will not be proved to be optimal, we will show that 
it rapidly outperforms the universal UI measurement. In 
addition we will show that in the case of continuous vari- 
able when the inputs are represented by coherent states 
the UI measurement can be also easily implemented by 
three beam splitters and two photo-detectors. 

Formally the unambiguous identification problem (UI) 
fits into the following framework: Three identical sub- 
systems A,B and C are prepared in unknown pure prod- 
uct states IV'?) A, and |V'i)s, |'02)c, respectively. Further- 
more, the subsystem A is guaranteed to be either in the 
same state as the subsystem B or as the subsystem C. 
Thus two types of states should be discriminated: 

|*1>ABC = WljA ® \'4'i)b ® |V'2)c , 

|«'2)aSC = |V'2)a® IV'OiJ «) 1^^2)0 • (1-1) 

The unambiguous identification machine is described by 
means of a positive operator value measure (POVM) 
consisting of three elements {iSg, ^'i, £'2}- Element Ei 
(respectively E2) corresponds to correct identification of 
l^'i) (respectively |^'2)) type of state and Eg corresponds 
to the inconclusive result. These elements must obey 
no-error conditions [Eq. p.2[ l] and constitute a proper 



POVM [Eq.dO])]: 

Tr[Eip2] = Tr[^2Pi] =0; p, = |*«)(*«| ; (1-2) 
E»>0,Eo>0; Eo + Ei+E2 = l. (1.3) 

It is assumed that the state (|^'2)) appears with a 
prior probability 7]i (772). The performance of the UI mea- 
surement is quantified by a probability of identification 
for a particular choice of states 

FdV-i), 1^-2)) = mTr[^iPi] + mTr[^2P2] ■ (1-4) 

Although PdV'i), IV'2)) is not a measurable quantity in 
the problems we consider, it will be very useful for com- 
parison of different UI measurements. Alternatively, we 
can use the average value to evaluate the performance of 
UI devices 

p{s) = 1 1 mi)M2)) dV'i#2. (1.5) 

In what follows we will denote the set of all pure states of 
a d-dimcnsional quantum system (qudit) by Sd and the 
subscript of P will indicate the used UI measurement. 
The optimality of UI measurements is defined with re- 
spect to their average performance, i.e. the aim is to 
optimize P{S). 



II. UNAMBIGUOUS IDENTIFICATION FOR 
QUBITS 

The solution to the unambiguous identification prob- 
lem for a qubit was given in Ref. [l3|. The optimal UI 
measurement depends only on prior probabilities 771 = 77, 
?72 = 1 — 77. Specifically, there are three different regions 
of values of 77 for which the POVM operators are specified 
as follows: 



0<77<i: 


Ef* 


= 0; 




E°P* 






Ef* 


- Alls 


® Vac; 


E°P* 




|<77<1 : 




= Is « 




E°P* 


= 0, 



where j^^^) = 1/V2(|01) - |10))Ai3, A = |(2-y^) and 

IV'^b) (V'asI ~ ''I'ab = ^^y^AB is a projector onto the 
antisymmetric part of the two qubit Hilbert space Tif^. 
The inconclusive result is associated with the POVM el- 
ement E°P* = 1 - E7* - E°P*. 



A. Relation to quantum state comparison 

Intuitively, the unambiguous state comparison and the 
unambiguous state identification are very closely related 
problems. Indeed, one can consider the following fam- 
ily of problems: given n + 1 qudits in unknown states 
-01, ... , ipn^ V'n-i-i- AH qudits [except the (71 -I- l)-th one] 
are guaranteed to be in different states. Decide whether 
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(n + l)-th qudit matches with one of the given qudits, 
or not. The unambiguous state comparison 0] is a task 
with n = 1 and unambiguous identification corresponds 
to n = 2 providing that the (n+l)-th qudit is promised to 
be in one of the states ipi, . . . jipn- Moreover, the identi- 
fication can be logically reduced to a series of state com- 
parisons, although such reduction does not necessarily 
give an optimal identification scheme. In particular, the 
UI measurement for n = 2 can be always designed out 
of the unambiguous comparison device in the following 
way: The experimentalist randomly chooses (with prob- 
abilities q and I — q, respectively) one of the reference 
states iV'c)) to be the input of the unambiguous 

comparator machine together with the unknown state 
\iPa)- Performing the state comparison Mab {Mac) 
the conclusive result means that the unknown state is 
iV'c) (IV'-b)) Formally this corresponds to a probabilistic 
switching between two measurement apparatuses result- 
ing in the UI measurement consisting of POVM elements 
= qlB^FfS^ E2 = {l-q)Ic®F%% Eo = I~Ei-E2, 
where F"^'^ denotes the POVM element associated with 
the conclusive result saying that two states are different 

0- . .. 

For a prior probability iji < i the optimal qubit 

UI measurement is projective, distinguishing the anti- 
symmetric and symmetric states of the subsystem AB. 
Moreover, the qubit C is not used at all, and the ex- 
change symmetry of the states with respect to systems 
A and B is measured distinguishing between the states 
\^i)a (8) \^i)b and \^2)a <8 \i'i)B- This is exactly the 
aim of a quantum state comparison of two unknown pure 
states. States of the type lip) (g) lip) are from the sym- 
metric subspace, therefore the projection onto the an- 
tisymmetric subspace unambiguously identifies whether 
the states of the subsystems are different. On the other 
hand projection onto the symmetric subspace is inconclu- 
sive, because both types of states \ip) (g) l^p), \ip) have 
nonzero overlap with it. Analogous considerations holds 
also for the interval r/i > | and a subsystem AC. For 
the aforementioned prior pro babilities the mean proba- 
bility of identification equals Popt{S2) = '7min/4, where 
'Hmin = min{77i, 772} G [0,1/5] and ^2 denotes the set of 
all pure states of the qubit. 

For equal prior probabilities the optimal measure- 
ment is a "true" POVM measurement, whose elements 
E°^*, Ej*'* are 2/3 times the above-mentioned quan- 
tum state comparison measurement elements ASym^f;^, 
ASym^^. In this case the mean probability of identi- 
fication Popt (5*2) is 1/6. Using the mixing strategy the 
success probability can reach at most 1/8. 



B. Unambiguous identification of equatorial qubits 

Let us consider a restricted set of pure states lying on 
the equator of the Bloch sphere \ip) = l/%/2(|0) + e"^\l)) 
with ip G [0, 2ti). Let us denote the subset of all equato- 
rial states by Seq- We are going to find the UI measure- 



ment, which optimizes the probability of identification 
Peq{Seq) avcragcd over the set Seq- Following the ap- 
proach used in Ref. [3 we obtain 



PeqiSeq) = mTr[El''ni] + 772rr [E^^fia] , (2.1) 
with average states 

= / / <<5j ® ® 'P2 d(pidip2 , (2.2) 

[^^) Jq Jo 

where we used the notation tpj = \Lpj){ipj\. After a little 
algebra this yields 

= ® {mm + |11)(11| + 2|7/;+)(V;+|)^B ; 

^2 = \^B ® (|00)(00| + |11)(11| + 2|V;+)(V^+|)^c , 

(2.3) 

with = -i=(|01) -I- |10)). We integrate the no-error 

conditions p.2p in the same way and obtain Tr[Ej'^r22] = 
rr[E2''r2i] ~ 0. Operators E^*, 17^ are positive therefore 
the previous equation means that the operators in the 
trace have orthogonal supports. Thus we determine the 
zero eigenvectors of the opeators fii and f22, f^2|ij) = 
and r2i|6j) = 0, respectively, where 

|ai) - |0)b ® IV-Ic) ; I&i) = |0)c®|^Ib); (2.4) 
W) = \\)b® I^Ic) ; 1^2) = |l)c ® IV'Ib) • 

These eigenvectors determine subspaces in which POVM 
elements Ej"^, E^« 

2 2 
Ef ^ ; = E PM{^^\ (2.5) 



can operate. Our goal is to maximize Peq{Seq), while 
keeping the POVM elements positive. Therefore we use 
equations (|2.3|) and (|2.5p to express equation (|2.1|) only 
in terms of the coefficients ajk, Pjk 



Peq{Seq) = ^{an + 022) + ^{Pii + P22) ■ (2.6) 



Accidentally at the same time the expression for Peq{Seq) 
coincides with P ppt ( 5*2) and the states \aj), \bk) are the 
same as in Ref. [13] [compare with Eqs. (3.19-3.22) of 
this reference]. Therefore the optimization task and the 
resulting measurement is in our case exactly the same as 
for the universal UI of qubits. As a result we see that 
the optimal UI measurement for equatorial states is the 
same as for the most general qubit state (specified at the 
beginning of this section). Hence, in this case the a priori 
knowledge does not help us to improve the performance 
of the UI measurement. 
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III. UNAMBIGUOUS IDENTIFICATION OF 
QUDITS 

A. The Swap-based approach 

The POVM elements for the optimal universal UI of 
qubits E^^'jEj^* are proportional to the projectors onto 
the antisymmetric subspace of the two-qubit subsystems 
AC and AB, respectively. The simple generalization of 
the aforementioned universal UI measurement to the case 
of qudits is the following POVM, which we abbreviate by 
sb (stands for the "swap-based"): 

Ef = cils ® ASym^c = ci^(l - SwapAc) ; 

Ef = C2lc ® ASym^B = ^-2^{^ - SwapAs) ; ^^'^^ 



E^ 



1 - Ef - Ef , 



where ASym^j^y denotes the projector to the antisym- 
metric subspace of X and Y particle. 

The positivity of Eg'', Ef'', Ef results in the conditions 
for ci, C2. Namely, from Ef'' > 0, Ef'' > we obtain that 
ci > and C2 > 0, whereas the inequality imposed by 
the positivity of Ef is not so apparent. Further we will 
calculate eigenvalues of Ef explicitly. Let {| be a 

basis of the qudit Hilbert space Hd- Then \ijk) = \i)A(^ 
\j}B ® \k)c is the basis of the three-qudit Hilbert space 
Ti.f^. The operator Eq'^ can be expressed in terms of the 
identity and the Swap operators so (ijfc|Ef |Zmn) = 
whenever {ijk} is not a permutation of {Imn}, i.e. ijk ^ 
a{lmn). In other words if we properly reorder the above 
basis the matrix Ef is block diagonal. The blocks arc of 
the following three types: 



the trivial 1x1 block (in I En 
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the 3x3 block with the matrix {ci (Mj)|Ef |(T2(iij)) 



• the 6x6 block 
((7i(ijfc)|Ef|(T2(uA:)). 



with the matrix 



The dimensionalities of the blocks are given by the num- 
ber of inequivalent permutations a of the three indexes. 

For qubits only the blocks of the first two types occur 
in the matrix Ef whereas for qudits (d > 2) blocks of 
all three types arise. Hence we reduced the problem of 
finding the eigenvalues of the positive operator Ef to an 
evaluation of the eigenvalues of the matrices mentioned 
above, which is treated in more details in Appendix A. 
It is shown there that the positivity of Ef imposes a 
particularly simple inequality 



Cl 



C2 < 1 



(3.2) 



Psb{\^l)M2)) = r;i(*i|Ef|*i)Ai3C+^2(«'2|Ef|*2)ABC 

•(i-i(vii^2)n. (3.3) 



regardless of the dimension d {d> 2) 
The probability (|1.4p reads 



Depending on the prior probabilities rji , 7^2 the values of 
ci,C2 maximizing Ps;,(|V'i)i l'02)) read: ci = 0, C2 = 1 for 
771 < 772 ! Cl = 1 . C2 = for rji > 7/2 and ci -|- C2 = 1 for 
Vi = ^2- Hence for equal prior probabilities the swap- 
based UI measurement is independent of the particular 
choice of ci and C2, i.e. 

Ps.(|^l),|^2)) = ^(l-|(^l|V^2)n. (3.4) 

However due to symmetry reasons we will further con- 
sider Cl = C2 = 1/2 in the case 771 =772, which gives 

Ef = ils® ASym^c, Ef = ^Ic® ASym^^ . (3.5) 

Averaging over the Bloch sphere and using the identity 
/ Isd K'/'il'/'2)pc''0ic^^2 = 1/d [in we obtain the average 
probability for the swap-based UI machines 



PsbiSd) 



I fd-l 



4 V d 



(3.6) 



Although the probabilities itself are independent of the 
dimension the average value converges to 1 /4 in the limit 
of d 00. This corresponds to an intuitive expectation 
that two randomly chosen unit vectors in Tid are more 
likely to be orthogonal for higher values of d. 



B. Optimal measurement 

Although POVM elements Ef'',Ef proportional to 
projectors onto the antisymmetric parts of the subsys- 
tem AC respectively AB, intuitively seem to be the best 
universal UI measurement, it was shown by A. Hayashi 
et al. [1^ that this is not the case for t^i = 772 = 77 = 1/2. 
The maximization of the average probability Popt{Sd) can 
be done by exploiting symmetry via joint representations 
of the unitary group U(E)U(E)U{UgU (d)) and the per- 
mutation group 5(3) permuting the subsystems A,B and 
C of ^ In our case Hayashi's optimal POVM mea- 
surement can be written explicitly in the form 



Ef* 
E°P* 



els ' 
etc ' 



' ASym^f, ; 



(3.7) 



where e — e\Tx is a specific operator. The parameter 
A specifies both U{d) and S{3) irreducible representation, 
Tx is the projector onto that invariant subspace of Tif^, 
and ex are non-negative real numbers. In our case only 
two irreducible U{d) representations specified by Young 
tableaux A = (2, 1, 0, . . . , 0), A = (1, 1, 1, . . . , 0) are rele- 
vant. The corresponding ca's are 2/3 and 1/2. Therefore 
we have 



U 



(2,1,0,. ..,0) 



(1.1,1 0)- 



(3.8) 



The projectors r(2,i.o,....o) a-nd 1 1 q) project onto 
the subspaces {Vs ffi Vas)'^ and Vas, respectively, where 
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Vs (Vas) is the totally symmetric (antisymmetric) sub- 
space of Hf^. Operators Ic ^ ASym^^, 1^(8) ASym^^ 
do not mix the subspaces {Vs © Vas)""" and Vas on which 
the operator e acts only a multiple of the identity, i.e. 

^l(Vs®v;4s)^ = 1. Therefore E°^* (analogously 



Ej^*) is essentially |lB(8)ASym^p except for Vas, where 
it is ^1_B ® ASym^^. Furthermore the POVM elements 
E°^*, Ej''* acting on the totally antisymmetric subspacc 
Vas do not contribute to Popti\ipi), \'ip2)) and Popt{Sd), 
because input states \'ij)ABC [ see Eq. (|l.ip ] are symmet- 
ric in a pair of subsystems. Thus for calculation of prob- 
abilities of identification we can as well use the operators 
El = f Is ® ASym^c = |Ef , E2 - f Ic ® ASym^^ = 



|E|^. Hence, in the optimal case 



PoptiHi),^) = ^(i-i(V'ii^2)n, (3.9) 



Popt{Sd) 



ld-1 
3 d 



(3.10) 



IV. 



UNAMBIGUOUS IDENTIFICATION OF 
COHERENT STATES 



Unlike previous sections, where we have considered 
an unambiguous identification of quantum states from 
a finite-dimensional Hilbcrt space Tid, here wc work with 
a semi-infinite dimensional Hilbcrt space of a linear har- 
monic oscillator Tioo, which models a single mode of an 
electromagnetic field (EM). The techniques presented in 
Sec. HI for qudits work for any dimension d. The result- 
ing POVM elements are expressed via constant multiples 
of projectors, which in large d limit define projectors on 
Ti'^. Therefore we have formally the same universal UI 
measurement also for states from Tioo and this measure- 
ment is optimal for the case of equal prior probability 

Our goal in this section is to show that an unambiguous 
identification of coherent states can be done with much 
better probability of identification than the optimal UI 
for all pure states from Tioo • The basic intuition for this 
is that coherent states form a very small subset Scoh of 
all pure states from Hoo and thus there could be a better 
way to identify them. The more reasonable motivation is 
based on the following observation. As it was mentioned 
in Sec. Ill Al for a special choice of the parameters 771,772 
the optimal POVM for qubits coincides with the optimal 
quantum state comparison measurement. Hence if there 
is a better quantum state comparison of coherent states, 
which can be used to design an UI setup for coherent 
states then this setup could perform better than the uni- 
versal UI measurement identifying all states from TYoo- 
E. Andersson et al. Q proposed such a quantum state 
comparison setup, which is also simply realizable by a 
beamsplitter and a photodetector. In what follows we 
will explain how their setup works. In addition, we will 
present a proof that it performs optimal quantum state 
comparison of coherent states. Then wc will show how 



it can be used to design an unambiguous identification 
setup for coherent states. 



A. Quantum comparison of coherent states 

A coherent state |a) is fully specified by a complex 
amplitude a. It is a pure state \a) £ Hoo, which is an 
eigenstate of the annihilation operator a\a) = a\a). A 
beamsplitter is a passive-optics device acting on a pair 
of EM field modes. Its action is described by the Hamil- 
tonian H = i6{ab^ — a^b) generating the unitary trans- 
formation U = e^'t{'^b''-a b) ^ where a,a^,6, 6^^ are creation 
and annihilation operators of the two modes. The op- 
eration of the beamsplitter is particularly simple for co- 
herent states and is determined by the interaction time, 
i.e. by the transmittivity T and the reflectivity R of the 
beamsplitter: 



\a) (g) 



Ta + y/R(3) 



Ra 



TP) (4.1) 



with T+R = 1. In comparison of coherent states we want 
to unambiguously distinguish between \a) = and 
I a) ^ 1/3). This is equivalent to distinguishing /3 — a = 
and /3 — a 7^ 0, which can be done by 50/50 beamsplit- 
ter (T = i? = 1/2) in the following way. The state of 
the second mode after passing through the beamsplitter 
will be either the vacuum |0) or the state ~ o;)) 

when a ^ p. Thus if we detect at least one photon 
in the second mode, which happens with probability 
1 - |(0|^(/3 - a))|2 = 1 - e-5l«-/5|', we are sure that 

the states were different. The detection of no photons is 
inconclusive, because all coherent states have a nonzero 
overlap with the vacuum. 

Now we are going to prove optimality of this setup serv- 
ing as unambiguous coherent state comparator. First we 
derive the optimal quantum state comparison measure- 
ment for coherent states. In general it is a POVM with 
two elements Hq''* (inconclusive result) and H^^* (unam- 
biguously indicating the inequality of states) obeying the 
following equations: 



Tr 



nria>x(a|®|«>y(a| 



= 



Va , 



H°f * > 0, Hf * > ; H°P* + Hf* = 1 







doTi 



(4.2) 
(4.3) 

over all coherent 
H7»(a| ® \a){a\ ^ Tr[H7*Ap.4) 



Integrating the no-error condition 
states we obtain the following condition 



which defines the operator A. The no-error conditions 
given by (|4.2p are completely equivalent to Eq. (|4.4p , be- 
cause the trace under the integral is nonnegative. The 
operators H^^* and A are positive and therefore Eq. (|4.4p 
means that these two operators have orthogonal sup- 
ports. Hence the largest possible support the operator 
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n°^* can have is the orthogonal complement to the sup- 



port of A. It is possible to show that 



(4.5) 



N=0 



where the vectors Ixjv) = 2"^ ElLo V (fe) \k)®\N-k) 
are mutually orthonormal, i.e. {xn\xm) = Sn^m The 
calculation of A is described in detail in Appendix B. 

Moreover the normalized operator is a projector 
with the same support as A. The support of the projector 
1 — - A is therefore the largest possible support of 11°^*. 
The optimal measurement must maximize the probabil- 
ity of revealing the difference of the states launched into 
the comparator: 



P«(|a),|/3))=Tr[nria)H®|/3)(/?|] 



(4.6) 



while keeping the positivity (0 < 11^''* < 1) and the no- 
error conditions satisfied. Combining these two condi- 
tions on the support of 11°^* yields H^^* = 1 — ^A. Thus 
the optimal unambiguous coherent state comparison is 
acomplished by the following projective measurement 



-A = V i\xN){xN\)xY ; 

TT ^ — ' 

= ^-Y.{\xn){xn\)xy- (4.7) 

Now it is sufficient to show that the described quan- 
tum state comparison setup performs the measurement 
()4.7|) . The mathematical description of the setup is sim- 
ple. First the beamsplitter acts on the modes X and Y 
via a unitary transformation U{9 ~ n/A) and then the 
photodetector discriminate between the zero-registered- 
photon result IIq = Ix ® (|0)(0|)y and the at-least-one- 
registered photon nf = 1x®{1- \0)(0\)y in the mode Y. 
Hence the measurement performed by the setup is given 
as 



(4.8) 



Unitarity of U guarantees that 11^* form a projective mea- 
surement as well as 11^. Therefore it suffice to show the 
equality Ilg* = IIq^*. It follows that 

oo 

= uH\N)x ® \0)yx{N\ ^ y{0\)U (4.9) 

N=0 

and the proof is concluded by showing that U^N)x ® 
\0)y = \xn)xy- The technical details are placed in Ap- 
pendix C. 

B. UI measurement with beamsplitters 

For the case of a qubit (Sec. II) we have seen explic- 
itly that the unambiguous identification is very closely 



related to the problem of the state comparison. In 
some sense the identification seems to consist of two 
state comparisons performed somehow simultaneously in 
a single run. However, because only a single copy of 
the unknown state is available, a very specific "cloning" 
machine should be used in order to make such reduc- 
tion of the identification problem possible. Usually the 
cloning machines |22| distributes the original quantum 
information (represented by quantum state) among sev- 
eral quantum systems Such approach results in a 
complicated entangled state such that the individual sys- 
tems are described by density matrices on average "clos- 
est" to the original quantum state. Unfortunately such 
cloning cannot be used for our purposes. The poten- 
tial clones described by mixed quantum states cannot 
be unambiguously compared with pure states. There- 
fore we need a very specific cloning machine producing 
the clones in pure states, i.e. l^p) (g) |0) |0)anciiia 
lip') (g) \tp') (g) |(/)(i/'))anciiia- The cloniug of coherent states 
has been analyzed in [2^, where it was shown that the 
single beamsplitter assisted by a linear amplifier is op- 
timal. Without the linear amplifier the beamsplitter 
alone performs on coherent states the transformation 
|a) g) |0) \VTa) g) I V^a), where R, T stands for refiec- 
tivity and transmittivity of the beamsplitter. And this is 
exactly a type of cloning we are looking for. Hence, the 
idea is to use one beamsplitter to clone the system A into 
two modes and afterwards use another two beamsplitters 
for particular state comparisons. Hence, in addition to 
the modes A,B,C we add an ancillary mode D set initially 
to vacuum, i.e. |$i„) = \a-?)A gi |q;i)b ^ \a2)c ® |0)d, 
where ja?) is guaranteed to be either \ai) or |a2). 

Our setup is composed of three beamsplitters the ac- 
tion of which is described by the unitary transformation 

1$™) ^ U3{DC)U2{BA)Ui{DA)\a,)\ai)\a2m (4.10) 

where Uj{XY) is associated with the j-th beamsplitter 
acting on the modes X and Y. The first beamsplitter Bi 
(with the transmittivity Ti) prepares two clones of the 
unknown state that is encoded in the mode A 



A- 



(4.11) 



The output system remains in a product state, hence the 
beamsplitters B2 and B3 can be analyzed separately 



B2 



\ai)B 



'T\oi-7 



^2^1 + /R^a7)^.12) 
- a/R2"i + \/T2Tia-?)A , 



and 

B-i ■■ I V^a?)D ® \a2)c 



I - V RsRia-! + \/T\a2)c- 



In case a? = ai we want that the beamsplitters B2,B3 
behave as in the comparison protocol of identical states 
|q;i), i.e. the modes A,C, respectively, should be 

transformed into vacuum. Such conditions tell us how 



the parameters of the beamsplitters should be adjusted, 
in particular, we obtain identities 



Ta = 



1 



1 + Ti 



2-ri ' 



(4.14) 



R, = 1. 



where we used the identity Tj 

The conditions specified by Eqs. (|4.14p can be met si- 
multaneously, therefore we set the transmittivities T2, T3 
accordingly. The final state of our four modes after pass- 
ing all three beamsplitters can be simply obtained from 
Eqs. (|4.12M.13P and reads 



\^out) = I V^(a? - ai))A ® IVT^ai + \/R2Tia9)B 
® 1^1^(02 -a7))c® I \/TVRia? + V^a2>r> • 

The field modes are still factorized and we can fo- 
cused only on states of modes A and C that are detecting 
whether the unknown state matches with ai, or a2. In- 
deed, depending on a? the modes A and C end up in the 
states 



a? 
a? 



(4.15) 



Measuring photon number in the modes A and C by 
photodetectors P2 and Pi, respectively, we can unam- 
biguously identify the unknown state. In each single run 
of the experiment we can distinguish four situations: i) 
none of the detectors click, ii) only Pi clicks, in) only 
P2 clicks, iv) both detectors click. In our situation both 
detectors cannot click at the same time, because at least 
one of the modes is in the vacuum. If only the detector 
Pi clicks from Eqs. (|4.15p we unambiguously conclude 
that a-? ~ ai. Similarly if only the detector P2 clicks we 
unambiguously conclude that a? — 0,2- If none of the 
detectors click we cannot determine which mode was not 
in the vacuum and therefore it is an inconclusive result. 

In the case a? = ai the probability of a correct identi- 
fication follows from equations (|4.15p and is given by the 
probability of detecting at least one photon in the mode 
C 



Pi = l-|(0|v^(a2-ai))|' = l 



(4.16) 

In case a? = a2 the probability of a correct identifica- 
tion is given by the probability of detecting at least one 
photon in the mode A 



^11 1 2 



P2^\- |(0|Vfi2(a2 - ai))r - 1 - 

(4.17) 

Thus the total probability of an identification for refer- 
ence states \oti) and \a2) is equal to 



Phs{\0il) M2)) = mPl + T]2P2 ^ 



(4.18) 



Next we want to optimize the performance of the setup 
by properly choosing the transmittivity Ti. The defi- 
nition of the uniform distribution on the set of coher- 
ent states is problematic, therefore we first focus on the 




FIG. 1: The beamsplitter setup designed for an unambiguous 
identification of coherent states. 



probability of identification for a particular choice of ref- 
erence states \oLiiB, \0L2ic expressed by Eq. (|4.18p . In 
fact, this will later help us to draw more general conclu- 
sions. By plotting the Pfcsdo^i), |q^2)) for various ranges 
of \ai — Q!2|,?7i G [0,1] and Ti S [0,1] one quickly finds 
that for the fixed values of 771 and |ai — aal the probabil- 
ity Pbsdai), \oi2)) is maximal for the values of Ti that de- 
pend on 771 and | ai — 02 1 . Thus in general for an arbitrary 
prior probabilities the optimal transmittivity Ti depends 
on the reference states to be identified. However we will 
show that in a special case of equal prior probabilities 
there is only one value of the transmittivity Ti , which is 
optimal for all reference states. This value turns out to 
be Ti = 1/2 as one would expect from symmetry argu- 
ments. In order to show this we calculate 

oi 1 

from Eq. (|4.18p for 771 = 772 = 1/2 and the condition for 
critical points (vanishing the first derivative) yields 



1 = 



(1 + Pi)^ 

(2-ri)2 



-|ai-0!2| (t 



-^1 



Tl ■ 
1 + Ti ' 



(4.19) 



For < Tl < 1/2 both terms on the right hand side (rhs) 
of Eq. (|4T9)) are greater than 1, for 1/2 < Ti ^ 1 both 
terms are less than 1 and for Ti = 1/2 both terms on the 
rhs are equal to unity. Thus Ti = 1/2 is the only critical 
point for all reference states and because of the second 
derivative being negative it is the global maximum of 
Pfcsdai), [02)) in the interval Ti g [0, 1]. 

Further, we will consider the UI of coherent states ap- 
pearing with equal prior probabilities. In this case the 
optimal choice of transmittivities for our three beamsplit- 
ter setup is Tl = 1/2, r2 = 2/3, T3 = 1/3. As a rcsuh 
we obtain that the probability of an unambiguous iden- 
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tification (|4.18p reads 

Pbs(|ai),|«2)) = l-e-^l"^""^l\ (4.20) 

C. Swap-like UI design for coherent states 

In this section we will propose different unambiguous 
identification measurement for coherent states. This ap- 
proach will be essentially the same as in the construction 
of the swap-based UI measurement for qudits, i.e. also 
motivated by the state-comparison problem. The dif- 
ference is that instead of considering all states we will 
be restricted to coherent states only, i.e. the role of 
an antisymmetric subspace is played by the projector 
^ ~ J2n=o \xn){xn\)ac, which was crucial for the op- 
timal state comparison of coherent states (|4.7p discussed 
in Sec. IV. A. Like before the conclusive POVM elements 
Ej^'^ ignore the identified mode and conclusively compare 
the states in the other two modes, i.e. reading that the 
states of these two modes are different. Such a POVM 
has the following structure 

oo 

Ef-^ = caB(^{t-Y.\XN){XN\)AC; 

N=0 
oo 

Ef-^' = C2ic® (1- 51 

El^f = l-Ef-^-Ef-^. (4.21) 

Let us fix the parameters ci, C2 and calculate 
Psbf{\cti),\ct2)) for this UI measurement. Using the 
identity J2n=o\xn){xn\ = f /c ItXtI «> |7)(7M7 [see 
Eq. (|4.5p ] we obtain 

Psbf{\ai),\a2)) = (7?ici +772C2) (4.22) 

X (1-- / |(ai|7)n(«2|7)Pd7). 

Taking into account that |(a|7)P = e^'" 'I and the 
rectangular identity |ap + |5p = |(|a + 6p + |a - 6p) 
we can express the integral / = /c |(cn|7)Pl(a2|7)prf7 
as follows 

/ = e~^l"i~"^l^ / e~^'"^+"^"^'^'^d7 
Jc 

^ J e-^h\'dj = |e-^l"i-"^l'(.4.23) 

Combining Eqs. (|4.22p and (|4.23p the unambiguous iden- 
tification probability reads 

P,b/(|ai), 1^2)) - imci + rnc2)il - e-5l«i-"2| = )(4 24) 

The positivity of the POVM elements Ef-'^, E^-^ is 
guaranteed by setting ci, C2 to be nonnegative. The op- 
erator E^'^ has a block diagonal structure in the basis 
consisting of number states ordered with respect to the 



increasing global photon number. Finding the eigenval- 
ues of this matrix is a difficult problem. Nevertheless, 
the POVM elements Ef-^ and £'2'''^ ai'c proportional to 
mutually overlapping projectors. Let \^j.)abc be a vector 
from supports of both these projectors. For example, 
can be a vector from the totally antisymmetric subspace. 
Then {^\Eq'''^ = 1 — ci — C2. Thus, the positivity 
implies that ci + C2 < 1. For equal prior probabilities 
'7i = '72 = 1/2 we obtain 

n6/(|ai),|a2))<i(l-e-^l— ^1^). (4.25) 

D. Comparison of UI measurements for coherent 
states 

In the previous sections we have discussed four differ- 
ent UI measurements that can be used to identify coher- 
ent states: 

i) the swap-based measurement, 

ii) the optimal measurement. 
Hi) the swap-like measurement, 
iv) the beamsplitter setup. 

The first two schemes unambiguously identify arbi- 
trary states of qudits in arbitrary dimensions. The re- 
maining two are designed to identify only coherent states. 
Although the comparison is usually understood in terms 
of average probabilities, we will adopt a different com- 
parison method evaluating the performance directly in 
terms of probabihties P{\ai), [0:2)) for all pairs of states. 
It turns out that for all the measurements these probabil- 
ities depend only on a scalar product of the states under 
consideration. 

As we mentioned at the beginning of this section qu- 
dit POVM elements E^*" and E"^* in the large-d limit 
define also POVM elements in Hoo- For simplicity we 
use the same notation for these operators. These two UI 
strategies are universal, so they work for any pure states 
from Hoc- If applied on coherent states the correspond- 
ing probabilities are given by Eqs. p.4p and (|3.9p 

P.bdai), 1^2)) = \{l |(ai|«2)|') = ^(1 - e-l"^-"^l') ; 

PoAl^i), |«2)) = i(l - I(ai|a2)n = ^(1 - e-l— ^1^) . 

In what follows we will compare Psbilc^i) j 
-Poptdai), |a2>), and Pbs{\ai),\<^2)), Psbf{\ai),\<^2)) 
which are probabilities of the identification in the UI 
strategies designed especially for coherent states [see Eqs. 
(|4.20p and ()4.25p ]. The following inequalities hold for ar- 
bitrary coherent states |ai) and |q!2) 

Psb <Psbf< Pbs ; Popt < Pbs ■ (4.26) 

Hence the same relations hold between the measurements 
also on average. They can all be derived in the same way. 
Let us define x = e-l"!-"^!" [x G [0,1]) and = ^. 
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FIG. 2: (Color online) The probability of identification 
P{\ai), \a2)) as a function of the scalar product (given by 
|ai — 02!) for four UI strategies applied on coherent states 
|ai), \a2). Starting from the bottom the two lowest lines cor- 
respond to universal UI measurements (the swap-based is in 
magenta (lowest line) and the optimal strategy is in black, 
respectively). The next two lines are associated with the 
UI measurements designed for coherent states (the swap-like 
measurement is in green while the and the three beamsplit- 
ters setup is in solid blue, respectively). The top (red) curve 
corresponds to the optimal discrimination probability among 
two known states. 



dance with the case iV = 2 we will firstly use a beamsplit- 
ter array to distribute (to "copy") the unknown state ja?) 
onto iV — 1 ancillary modes (which are initially set to the 
vacuum). After this redistribution of the information we 
will simultaneously implement A^-fold state comparison 
to unambiguously identify the unknown state. 

The quantum state distribution can be done with TV — 1 
beamsplitters Bj (described by parameters Rj, Tj) acting 
on the j-th ancillary mode and the mode of an unknown 
state. The bcasmplitters are applied sequentially split- 
ting the unknown state into N modes (see Fig. ^ so that 
each of them end up in the state | -^a?) . After a little al- 
gebra one can derive the following values for reflectivities 
and transmittivities of the j-th beamsplitter 



T, 



N-j 
N-j + 1 
1 

N-j + 1 



(4.27) 



Altogether these beamsplitters will implement the trans- 
formation 



|0) 



(giAf-l 



1 



N 



(4.28) 



All the probabilities are zero for ai = a2 {x = 1) as 
they should, because in such case reference states coin- 
cide. The validity of these inequalities can be proved by 
showing the reversed inequalities for the first derivatives 
of probabilities (|4.26p with respect to x, i.e. 



After this transformation is completed we will simul- 
taneously apply N beamsplitters performing the quan- 
tum state comparison of states and \aj). Let us 
denote by Cj the beamsplitter comparing the unknown 
state with the state of the j-th mode. Each of them 
performs the following transformation 



dxPsb > dxPsbf > dxPbs ; dxPopt > d^Pb. 



> 



> 



> 



3x3 



3x3 



The second row of inequalities obviously holds in the in- 
terval X € [0,1], so the inequahties (|4.26p are proved. 
More quantitative insight is given in Fig. ([2]) showing the 
dependence of the probability of identification for the 
considered UI measurements on the value of \ai ^ a2\- 

As a result we can conclude that the beamsplitter 
setup designed for an unambiguous identification of co- 
herent states performs better than other devices includ- 
ing the optimal universal UI measurement. Another re- 
markable feature is that the beamsplitter setup attains 
^'fc.sdo^i), \ct2)) ~ 1 for large values of \ai — 02], i.e. in 
the limit when two coherent states are orthogonal. 



E. Unambiguous identification of N reference 
states: A quantum database 

Now let us consider the problem of unambiguously 
identification among N coherent states |ai), . . . , |aAr). 
Our aim is to modify the proposed beamsplitter scheme 
to address this slightly more general problem. In accor- 




where we used = ■j^,T;: = 
\a, (3) tor \a)(E)\(3). As before the photodetcctors monitor 
the photon number only in the modes originally in states 



and the notation 



\aj) that at the output arc in the state 



1 



(a?-«j))]<»l 



1 



(aAT-a?)) . (4.29) 



Hence if ja?) = \ak} the fc-th mode is in the vacuum 
state and all other modes are excited (i.e., populated by 
photons). Therefore if all photodetcctors except a single 
one click then we can unambiguously conclude that the 
unknown state matches with the initial state of the mode 
corresponding to the detector which did not fire. For all 
other combination of outcomes the result is inconclusive. 
This implies that the probability of unambiguous iden- 
tification for the reference states |q!i), . . . , |aAr) is equal 
to 



Pi 



\ai) 



N 

E 

.7=1 



N 

^.11(1 



=j\ak- 



01.30) 
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FIG. 3: Unambiguous identification measurement setup iden- 
tifying among A*" colierent states. 



If we set j]j = (for all j) and ak — ae^^'^ then the 
probability can be simplified to 



P(|ai),...,|a^))= n 



fc=i 



= (l_cos(^fc)) 



(4.31) 



Concluding this section we note that the unambiguous 
discrimination among N unknown states described above 
can be considered as a search in a quantum database com- 
posed of N elements, i.e. N different though unknown 
coherent states \aj) that are encoded into N modes of 
an electromagnetic field. We point out that we have only 
a single copy of each of the states \aj) so one can not 
acquire a complete classical knowledge about the state. 
This set of N states corresponds to a quantum database. 
In addition we have the (TV -|- l)-st mode of the light 
field in the state |q;?). The search of the database corre- 
sponds to the task of matching of two modes such that 
ja?) = \aj). So we can say that the two modes are in the 
same state without knowing what the state actually is. 



V. CONCLUSION 

In this paper we have addressed the problem of an 
unambiguous identification of unknown coherent states. 
We have explicitly designed two UI measurements tak- 
ing into account the a priori knowledge about a partic- 
ular family of states and compare these measurements 
with the universal unambiguous identification, i.e. the 
UI measurements (either the swap-based or the optimal 



one) that can be applied for all states. Our main goal was 
to design a simple experimental setup consisting of three 
beamsplitters (see Fig.l) that performs best. Finally, 
we have generalized the problem and proposed the mea- 
surement unambiguously identifying among N coherent 
state. 

The beamsplitter setup was motivated by an intuitive 
reduction of the unambiguous identification problem into 
specific "cloning" task and an unambiguous state com- 
parison. We have proved that the state comparison orig- 
inally proposed in is indeed the optimal one (this was 
implicitly conjectured in Ref . 0] but was not proved) . It 
is interesting to compare UI measurements described in 
this paper with an UI measurement given as a mixture 
of two optimal unambiguous comparison measurements, 
i.e. El = qlB®Ff^, E2 = qIc®F%, Eo = I-E1-E2. 
Let us consider the UI problem for pure states belonging 
to a set S such that \ip){Tp\ (g) {i>\dip = Il{S) is a 
projector. Under such assumption the conclusive result 
for an optimal unambiguous comparison of states from S 
is associated with the positive operator F'^^^ = I — 11(5). 
li S — Sd is the set of all pure states of d-dimensional 
system then F'^'^ = / - Il{Sd) = I - Hsym = ASym. 
Hence, the corresponding POVM consists of operators 
El = \Ib ® ASym^c', E2 = \lc ® ASym^^, i.e. (for 
d > 2) it is exactly the SWAP-bascd UI measurement 
as specified in Ea. (|3.5p . Similarly for coherent states we 
have F'^'^ = I - n(S'coh) = / - |A (see Eq.gl])) and 
the mixing of optimal unambiguous comparison strate- 
gies results in the same POVM as in Eq. (|4.2ip . Based on 
our results we can say that although for purposes of the 
state comparison problem the operator 1 — 11(5) is the 
optimal solution, the resulting UI measurement based on 
mixing of such optimal unambiguous comparison strate- 
gies is not the optimal one. 

The proposed beamsplitters setup (see Fig.l) for un- 
ambiguous state identification can be compared with the 
measurement proposed in Ref. Q discriminating opti- 
mally among two known coherent states. Both of them 
consists of three beamsplitters, but arranged differently. 
An interesting observation is that the differences between 
the probabilities are not very large (see Fig. 2) and even 
more surprising is the fact that two unknown nearly or- 
thogonal coherent states can be identified almost per- 
fectly. For the universal optimal UI measurement (see 
Fig. 2) there is a significant gap between the probabilities 
for state discrimination and state identification. 

The proposal of unambiguous identification of coherent 
states is extended for an arbitrary number N of reference 
states. We proposed an experimental setup consisting 
of 2N — 1 beamsplitters for unambiguous identification 
among N unknown coherent states. This setup can be 
considered as a search in a quantum database. The el- 
ements of the database are unknown coherent states en- 
coded in different modes of an electromagnetic field. The 
task is to specify the two modes that are excited in the 
same, though unknown, coherent state. The analysis of 
this aspect of unambiguous identification problem is be- 
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yond the scope of this paper and details will be presented 
elsewhere. 
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For qudits, d > 2, at least one 6x6 block appears in 
the matrix of Eg**. The eigenvalues of the 6x6 block 
satisfy the following inequality: 

\(6) ^ 1,(6) ^ i(6) ^ <(6) 

so the sufficient condition for positivity of Eq^ when d> 2 

is A^^^ = 1 - Ci - C2 > 0. 

APPENDIX B: CALCULATION OF A 



APPENDIX A: EIGENVALUES OF E'o 



In this appendix we will explicitly evaluate the opera- 



tor 



The operator Eg'' is defined in Eq. (|3.1|l . As we have 
already mentioned in Section IIIA this operator is block 
diagonal and consists of three types of blocks. 

1. Trivial {iii\W^\m) = 1. 

2. 3 X 3 matrix {ai{iij)\W^\a2{iij)): 



iij l-ci/2 ci/2 
iji 1 - C2/2 C2/2 
Kjii ci/2 C2/2 l-ci/2-C2/2y 



with eigenvalues 



A 



(3) 



1: 



, (3) _ 2 - Ci - C2 ± \/cl- C1C2 + c 

-^2,3 - ;s 



3. 6 X 6 matrix (CTi(ijfc)|EgV2(ijfc)>: 



ijk X 

kji ci/2 

jki 

ikj 

kij 

\jik C2/2 



kji 
ci/2 

X 
C2/2 









jki 


C2/2 

X 
ci/2 







ikj 



ci/2 

X 
C2/2 





kij 




C2/2 

X 
ci/2 



jik \ 
C2/2 







ci/2 
X ) 



(Al) 



(A2) 



(A3) 



\o)x®\oi)Yx{iA®Y{pL\ da. (Bl) 



First expand the coherent state \a) in the number basis 
and use the polar paramctrization of the complex plane 



00 I. 



a 



2 —^\k), a = re"*, da ~ rdrdcj) 
in order to write 



-2| 



.1^ ^ 



E 



k,l.m,n—0 

1 

\/kU\m\nl Jo 



\/kW.mln\ 



\k){m\ (g) \l){n\da 



^-2r^ ^k+l+m+n+l 



dr 



.k, 1,771,71 — 

X / e''^'^''+^-"'-"U(j)\kl)(mn\ 
Jo 



The integration over (/) yields 2T:dk+i,7n+n and the inte- 
gration over r after a simple substitution and usage of 
Sk+i,7n+7i leads to r(fc + I + l)/2''+'. Using the identity 
r(A: + ; + 1) = ^(fc + /)!^(m + 77,)! we obtain 



A 



kJ.m.n—O 



Jk-i-l,77l-\-7l 



\kl) {mn\ 



where X = 1 — ^ ^ The corresponding eigenvalues 
read 



a: 



(6) 

v(6) 
^5,6 



1 
1 

2 — Ci— C2 + -^ c'j — CiC2+cl 



(A4) 



2 — Ci— C2 — -^cj — CiC2+c| 



For qubits we have 2'^-dimensional Hilbert space and 
Eg* is represented by 8 x 8 matrix with two 3x3 blocks 
and two 1x1 blocks. So for qubits the sufficient condition 
for positivity of Eg'' reads 



Ci - C2 ± \ C 



C1C2 



Defining a new summation index N = m + n = k + l and 
using Sk+i.rn+n wc eliminate one sum to obtain 



N 



A 



\ k, / 

if 2* 



00 
N=0 

where we defined mutually orthogonal vectors 

N 



\k)\N -k){m\{N -m\ 



(B2) 



k=0 

i-e. {xn\xm) = Sn,m- 



N 



\k)x®\N~k)Y, (B3) 
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APPENDIX C: CALCULATION OF U^N)x ® I0)y state \xn} = U^N)x\0)y 

The action of the 50 / 50 beamsphtter unitary transfer- ^ 
mation U on the creation operators , b'' of modes X \xn) = (U^a^U)^ 



and Y is given by 



1 



U^a^U = (CI) 
v2 



(at+fe^)^|0)| 



'iV! 



1_ A lk\{N-ky.) fN\ (at)'^ (&t) 
V2^ " - . zf 



(7tfet[/ = i=(_at+6t) (C2) 2f^„V M \k J ^{N - k)\ 



N 



Let us express the number (Fock) state via the ere- — "^^^ J [^\k)\N ~ k) . 

ation operators acting on the vacuum, i.e. \N)x = 2"2" y V^/ 

(a^)^/%/M|0)x- The vacuum is unaffected by both U 
and W . Having aU this in mind we can calculate the 
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